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ABSTRACT 

The X-ray holes at the centre of the Perseus Cluster of galaxies are not all at the same po¬ 
sition angle with respect to the centre of the cluster. This configuration would result if the 
jet inflating the bubbles is precessing, or moving around, and the bubbles detach at different 
times. The orientations which best fit the observed travel directions are an inclination of the 
precession axis to the line of sight of 120° and an opening angle of 50°. From the timescales 
for the bubbles seen in the cluster, the precession timescale, Tprec, is around 3.3 x 10^ yr. The 
bubbles rising up through different parts of the cluster may have interacted with the central 
cool gas, forming the whorl of cool gas observed in the temperature structure of the cluster. 
The dynamics of bubbles rising in fluids is discussed. The conditions present in the cluster are 
such that oscillatory motion, observed for bubbles rising in fluids on Earth, should take place. 
However the timescale for this motion is longer than that taken for the bubbles to evolve into 
spherical cap bubbles, which do not undergo a path instability, so such motion is not expected 
to occur. 
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1 INTRODUCTION 

High resolution radio observations of a number of extragalactic ra¬ 
dio sources show that the opposed jets contain inversion symme¬ 
tries, which pr ovides evidence th at the central e ngine is precessing , 
e.g. NGC 326 iEkers et alJl978^ and NGC 315 iSridle et alJ 19791) . 
In a review of quasars and radio galaxies which exhibi t bent jets, of 
which some are inversion-symmetric about the core. iGower et alJ 
^1982^ present ten sources whose morphology can be explained us¬ 
ing a precessing jet with periods ranging b etween 10^ —> 10^ yr. 
The e xtraordinary galactic object SS433 iHiellming & .TohnstonI 
19^^H succe ssfully described with a precessing jet model by 
Margonl il984|), and has been the subject of a recent study by 
Bfunddl & Bowl^ FOOSl) . 

The radio source at the centre of the Perseus Cluster, 3C84, 
is an extended source and shows an “S” or “Z”-shaped morphol¬ 
ogy at GHz frequencies. These lobes have been observed to anti- 
correlat e spectacularly with the observed X-ray em ission from the 
cluster iBdhringer et alJll9^ Ipabian et alJl200(l) . These depres¬ 
sions in the X-ray emission are thought to be the result of the 
jet from the central engine blowing bubbles of relativistic plasma 
into the thermal Intra-Cluster Medium, pushi ng the ICM aside 
iGull & NorthoveJl97.3l : lBdhringer et alll993h . When their buoy¬ 
ancy velocity is greater than their expansion velocity, the bubbles 
are expected to detac h and rise up buoyantly through the ICM 
iChurazov et alJEoOd) . This is observed in the Perseus Cluster, 
to the North-West and to the South of the core where there are 
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depressions in the X-ray emission which do not have any asso¬ 
ciated GHz radio emission; however, at lower frequencies (e.g. 
330 MHz) “spurs” of emis sion are observed po inting towards these 
so-called “Ghost bubbles” iPabian et ai]200a). 

The two pairs of bubbles which are clearly visible in the X-ray 
emission from the cluster belong to two different “generations,” the 
young, radio-active pair and the older “ghost” pair. The bubbles are 
observed at different position angles with respect to the radio core 
in the centre of the cluster. In Section|^we discuss the morphology 
of the radio emission at the centre of the cluster, linking in other 
features in Section^ Section|3outlines the model for the precess¬ 
ing jet, the resultant parameters of which are presented in Section 
13 We calculate the timescales for the precession in Section|3and 
discuss the model in Section|7| We put forward two candidates for 
causing the precession in Section [3 and discuss the dynamics of 
rising bubbles in Section|3 


2 THE RADIO BUBBLES OF THE PERSEUS CLUSTER 

The extended radio emission from the radio source does not de¬ 
crease smoothly from the core; th ere is a sharp drop in radio flux 
in the southern lobe (see Fig. m& IPedlar et a 011993). This could 
be because there are two bubbles present, an older one which is in 
the process of rising buoyantly and a younger one which is still at¬ 
tached to the jet and being inflated. The X-ray emission also shows 
features in the region where the drop in radio flux occurs, there is 
a swath of excess X-ray emission which could arise from the ICM 
flowing in behind the recently detached bubble (Fig.0. The shape 
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Figure 1. LEFT: The X-ray emission from the centre of the Perseus Cluster from the 200 ks obsei'vation of iFabian et alJl^03s^ ■ corrected and accumulatively 
smoothed with a ct = 15. RIGHT: A 330 MHz radio map of 3C84. The shapes used for the different hubbies and the lines joining their centres to the radio 
nucleus are shown along with their names. 


of the X-ray decrement corresponding with this “Southern Outer” 
bubble is also similar to that of a spherical cap bubble. 

There is no equivalent sharp drop in the radio flux in the north¬ 
ern lobe, which, if there are also two bubbles present in this lobe, 
may result if younger one is behind the older one. This could soften 
the change in flux between the two bubbles as projected on the sky. 
In the case that the northern and southern radio lobes contain two 
bubbles each, we identify them as shown in Fig.Q A spectral index 
map of the central regions of Perseus (Fig.|3 shows that the bub¬ 
bles we have labelled as Outer have a steeper index than the Inner 
ones (—1.7 vs. ~ —1.2 for Sv oc iz“). This supports the idea that 
these Outer bubbles are older than the Inner ones, as the electrons 
have aged, steepening the spectrum. 

What can be seen from these observations is that the bubbles 
do not necessarily travel in the same directions when they detach. 
This is most obvious from the ghost bubbles. This may be due to the 
bubbles being somehow deflected from their initial path. Another 
explanation may be that the jet produced by the Super-Massive 
Black Hole (SMBH) is precessing, and the bubbles detach at differ¬ 
ent phases in the precession cycle. This analysis assumes that the 
bubbles travel outward from the centre in the direction the jet was 
pointing in at the time that they detached, with the precession time 
longer than the bubble inflation time. The current positions of the 
radio bubbles allow some constraints to be placed on the orientation 
and timescales of the precession. 

However, Leonardo da Vinci ob served that bub bles rising in 
liquids can move in a zig-zag fashion La Vincill.51.5h . and this may 
also be an explanation for the difference in observed position an¬ 
gles of the bubbles. However the analysis in Sectionl^leads to the 
conclusion that this behaviour is not relevant in the Perseus Cluster. 

The northern ghost bu bble has been likened to a spherical cap 
bubble rising through fluid iFabian et all20033l . The shape of the 
bubble is very similar to such a bubble in profile, and so can be rea¬ 
sonably expected to be travelling out from the centre of the cluster 
in the plane of the sky. Observations of the Hq filaments which are 
found around NGC1275 match some of the flow lines expected be- 
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Figure 2. The 330 to 1500 MHz spectral index map of the Perseus Cluster 
(G. Taylor, private comm.) showing the shapes for the bubbles. 

hind such a bubble very closely. iHatch et*^ ] <20Q5t show that the 
line of sight velocities for the Horseshoe-filament which sits behind 
the North-Western Ghost bubble are very small for the sections of 
the filament that point towards the centre of the cluster, implying 
that they probably move at ~ 90° to the line of sight. 

Very Long Baseline Interferometry (VLBI) maps of the region 
close to the core show an asymmetry between the southern and 
northern jets. While the southern jet is seen at all radio frequencies. 
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the inner parts of the northern jet are not detected below 15 GHz. 
At the lower frequencies free-free ab sorption has been s hown to be 
a reason why not all of the jet is seen iWalker et al J200nll . At higher 
frequencies, the inclination of the jet would play a role, with the 
southern one being inclined towards us. The amount of beaming, 
and hence asymmetry depends on the jet’s bulk velocity. This in¬ 
clination of the jet matches the orientation of the larger scale lobes; 
the northern Inner on may be behind the older Outer one. 

Measurements of t he VLBI jets can give so me constraints on 
the current jet direction. iKrichbaum et aljil992h calculate that the 
jet is at $1 14.4° to the line of sight, and even $1 2.7° close to 
the core for the milliarcsecond jet; however they find the jet to be at 
39.4° ^ a ^ 58.2° to the line of sight when observed at arcsecond 
scales, requiring it to bend by 25 — 45°. We take the jet to be at 
~ 10° to the line of sight. 


3 OTHER FEATURES 

3.1 Low Temperature Whorl and Ancient Bubbles 

A spiral struct u re in the Perseus Cluster was first noticed by 
Qiurazo^^jj_j2000|) and is seen in the tempe rature maps of 
Sande^Tt alJ i20oJi iiFig. 1^. Ipabian et alJ JTQOcIi suggested that 
this whorl could arise from the merger of a cool subcluster, as strip¬ 
ping of gas as it passed through the ICM would have left a wake. It 
could also form as the result of cooling gas falling into the centre of 
the cluster if the gas started with some initial angular momentum. 
However, bulk motions of the gas would be expected to have an 
effect on the Hq filaments ISection f3.2> and as they are linear this 
is unlikely as an explanation. If the jet is precessing, then the bub¬ 
bles rise through the ICM at different position angles. They would 
push low temperature material in front of them or drag it up in their 
wake as they rise. The North-Western Ghost bubble falls on the in¬ 
ner edge of this whorl. There is also a hint that there is an even older 
(“ancient”) bubble to the North, seen as a crescent-shaped bump on 
the outer surface of this whorl (3hl9m45s, 41:33°) (Fig.|3. If this 
is an ancient bubble, then the decrement in the X-ray emission is 
now so small that it cannot be clearly seen in the X-ray images 
of the cluster. However the change in the temperature of the ICM 
caused by its motion up through the cluster is still visible. 

This whorl is most easily explained with a precession model 
where the precession axis is face-on to us with a large 
90°) opening angle. As the bubbles detach from such a slowly 
clockwise-precessing jet they push up cooled, low temperature gas 
from the centre of the cluster. The next bubble lifts material from a 
slightly different position angle. This low temperature gas could be 
connected together by magnetic fields present in the cluster. Thus 
the portions of gas which are not lifted/pushed up by the bubbles 
are pulled up by the portions of gas that are, and so a whorl forms 
as opposed to discrete portions of gas around the rising bubbles. 

The interaction of the Northern bubbles (Ancient and Ghost) 
with the whorl is fairly clear, the Southern Outer bubble also sits 
behind it and so could just be starting to lift this section of the 
low-temperature gas out of the centre of the cluster. There is some 
extension of this cool gas to the West in the same dire ction as a 
high abundance ridge presented in ISanders et al.l ^OO^). This has 
also been interpreted as the remnants of an ancient bubble which 
lifted up the enriched material from the centre of the cluster. The 
interaction of the Southern Ghost bubble with the low temperature 
gas is not as clear, but it sits just South the coolest gas at the centre 
of the cluster. It may have lifted cooler gas out from the centre. 
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Figure 3. A smoothed version of th e tem perature map 
of the Perseus Cluster iSanders ~ aD |200^ . interpolated 
with the NATGRID natural neighbour interpolation library 
I http://ngwww.uc ar.edu/ngdoc/ng/ngma th/natgridynnhome.html i. For 
more details see iFabian et alJ llOOSJl . There are some edge effects in 
the image resulti ng from the smoothin g. White circles are excluded point 
sources, listed in ISanders et alJ ^2QQ4^ . The approximate locations of the 
older bubbles present in the cluster are shown. 

but if it were travelling close to the line of sight, in projection the 
distance the gas would have appeared to have been moved is not 
that far. 

It would be expected that if this whorl is the result of the 
jet/bubbles interacting with the ICM, then two spirals should be 
seen, one resulting from each jet. However, it seems as if there is 
only one. The explanation for this could be that the cooling has 
resulted in an asymmetric arrangement of gas. The northern bub¬ 
bles could then have broken through this gas, whereas the southern 
ones have detached in such an orientation that the spiral appears 
continuous. 

The whorl may, of course, be the result of a chance projection 
of two rings; a larger outer one which forms the northern and south¬ 
western parts of the whorl and an inner one around the northern 
outer bubble and the eastern part of the southern younger bubbles. 
We do not offer an explanation of the formation of this structure 
but merely put forward than the the whorl may not be one single 
structure. 


3.2 Ha Filaments 

The large Ha nebul a around NGC 1275 has been st udied at a vari¬ 
ety of wavelengths (MinkowskilllQ.STtlLvndsUlQTch . and deep im¬ 
ages obtained by the WIYN telescope show ext ended filaments 
which stretch over 50 kpc fr om the central galaxy (Conselice et alJ 
l200ltl . iFabian et alJ (20031^1 showed that there was a correspon¬ 
dence between the soft X-ray and the Ha filaments. They also 
showed that the filaments may represent the flow lines of the ICM 
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Table 1. Bubble directions and Timescales. 

Angle from North (0°)^ Timescalest/10^ yr 
Lobe NS NS 


Jet 

0 

180 

- 

- 

Inner 

15 

205 

1.1 

1.0 

Outer 

355 

145 

1.5 

1.9 

Ghost 

55 

190 

7.5 

7.1 

Ancient 

15 

130 

10.0 

10.0 


t The angles are measured Clockwise when looking at the figures pre¬ 
sented here. 

i The timescales (ages) have not been adjusted for projection effects. 
Note: All angles and timescales are approximate. The Ancient bubbles are 
discussed in Section|^3 

in the wake of passing bubbles, for example the “horseshoe” fila¬ 
ment which sits behind the North-Western ghost bubble. There is 
a large collection of filaments stretching up northwards from the 
cluster in the direction of the possible ancient bubble mentioned 
above. If the northern ancient bubble is a real feature and the fil¬ 
aments are linked with it, then they are likely to ha ve been drawn 
up out of the cluster by the passage of the bubble. ISanders et alJ 
| 200 ^ also find evidence for filaments pointing in the direction of 
the South-Western ancient bubble. 


3.3 Spiral Arms in NGC1275 

In Hubble Space Telescope (HST) images of the nucleus of NGC 
1275, ripples were noticed in the galaxy light iHoltzman et alJ 
119921: ICarlson et ID 1 19981) which resemble the arms of spi¬ 
ral galaxy. These ripples spiral in the same sense as the low- 
temperature whorl seen in the X-ray gas (Fig. |3j and so may be 
connected. If there is a processing black hole at the centre of the 
galaxy, then the precession of the jet may be linked with these fea¬ 
tures. 


4 PRECESSING JET MODEL 


From the shape of the low temperature whorl, the precessing jet 
may arranged such that the precession axis points close to the line 
of sight, with an opening angle of around 90° (wide cone), or, as 
there is some reflection symmetry about the East-West axis a nar¬ 
row cone precession where the precession axis is aligned North- 
South with a smaller opening angle. 

To describe the precessing jet, the model outlined in 
iBarker & Bv73 ll98ll) has been adapted and is shown in Fig .|3 
where 7 is the angle the precession axis makes with the line of 
sight, and S is the angle the jet makes with the precession axis. The 
angle the jet makes with the line of sight, a, is given by: 

cos a = cos 7 cos (5 + sin 7 sin (5 cos ( 1 ) 


where (f> is the phase angle of a clockwise rotating jet, and is 
0° when the jet points towards Earth. When a = 90° then 
cos (j> = — cot 7 cot 5. 

If tp is the angle between projected direction of the jet on the 
plane of the sky and the projected direction of the precession axis, 
then 


, sin (5 sin A 

^-X--1- 

Sin 7 cos d — cos 7 sin d cos (p 


( 2 ) 


An extreme value of tan pi occurs at cos (j) = tan 5 cot 7 . 

The angles from North (clockwise as seen in the Figures) and 




PLANE OF THE SKY 


Figure 4. The model of the precessing jets. 


the timescales for all the bubbles are listed in Table Q We have 
used the sound speed timescale for the Inner and Outer bubbles 
and the buoyancy timescale for the Ghost and Ancient bubbles. 
The buoyancy timescales estimate the time taken for a bubble of 
the currently observed size to rise buoyantly from the centre of 
the cluster through a uniform medium of that it is currently ob¬ 
served to be in. The sound speed timescale estimates the time 
taken for a bubble of the currently observed size to be blown if 
the expansion is subsonic. These two timescale measures are dif¬ 
ferent, but provide the best estimate of the age of the bubbles at 
different stag es in their formation (for further discussio n of these 
timescales see lDunn~ Fa^l 2 ^;|D unn et a 011^). The dif¬ 
ference in timescales between the Inner and Outer bubbles is small 
(0.4 —> 0.9 X 10^ yr), whereas the differences are larger between 
the Outer/Ghost and Ghost/Ancient bubbles. This may be a result 
of the different estimates of the timescales that have been used. 
The buoyancy timescale of a bubble, being a rise time, is gener¬ 
ally larger than the sound speed timescale, which is an expansion 
time, and so the jump in ages between the Outer and Ghost bubbles 
may be a result of this. It is unlikely that the differences are largly 
the result of anisotropies in the ICM, as th e pressure exerted o n 
the bubbles changes only by a factor of ~ 2 iSanders et all200'3) . 
Projection effects may remove (or enhance) this jump as the sound 
speed timescale is unaffected by the orientation, being an expan¬ 
sion timescale; whereas the buoyancy timescale does depend on 
the inclination of the rise-path of the bubble to the line of sight. 


5 CALCULATION OF THE PRECESSION PARAMETERS 

When calculating the best fitting parameters, the angle of the jet to 
the line of sight has not been used as a constraint as the actual angle 
is uncertain. However, as the difference in the emission from the 
Northern and Southern parts implies that the southern jet is coming 
towards us, this has to be satisfied by the parameters. 

As the Northern Ghost bubble appears to be moving close to 
the plane of the sky we assume that a = 90° for this bubble. The 
shape of the Southern Outer bubble also appears to be such that it 
is also moving close to the plane of the sky; this has also been used 
to constrain the possible angles. However, as there are only two 
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Table 2. Phase angles and adjusted timescales 



Phase angle,^ 0 

Angle to 1.0.s. 

Timescales^ 

Lobe 

N 

S 

N 

S 

N 

S 

Jet 

182 

182 

170 

10 

— 

— 

Inner 

6/179 

19/177 

70/170 

108/10 

1.1 

1.0 

Outer 

341 

293 

72 

94 

1.5 

1.9 

Ghost 

67 

0 

94 

110 

7.6 

8.1 

Ancient 

6 

247 

70 

55 

11.4 

15.3 


t The phase angle is measured for the Northern jet, and the precession is 
assumed to be clockwise as see in the images presented here, 
i The timescales have now been adjusted for projection effects. They are 
accurate to around 20 per cent. The sound speed expansions timescale has 
been used for the Inner and Outer bubbles, the buoyancy rise time for the 
Ghost and Ancient bubbles. 

Note: All angles and timescales are approximate. Timescales in 10^ yr. 

bubbles which provide constraints on a, and these appear to move 
in the plane; a close to face-on precession axis with a high value 
for 5 (around 90°) would also fit the data. 

Although this arrangement is a possiblity, such an extreme 
opening angle would be very unusual for a black hole of this mass. 
Also it would be expected that the bubbles would have detached at 
all position angles over the course of a few precession timescales. 
However the position angles of the northern bubbles are very simi¬ 
lar to those of the southern bubbles, almost mirror symmetry about 
an East-West axis. 

From EquationsQandl^and the observed projected angles on 
the sky, the best fit occurs when the projected angle of the pre¬ 
cession axis on the sky is 10° (clockwise of North as seen in Fig. 
0. Then there are various combinations of 7 and 5 which fit the 
observed data. However, we choose those with the smallest S (pre¬ 
cession opening angle) are 7 = 60° or 120° with S = 50° (±5°). 

To decide which of the two scenarios is the more plausible we 
look at the VLBI jet. The Southern jet is brighter than the Northern 
jet, which implies that it is Doppler beamed, and so is thought to 
be travelling towards us. The angles of the jet to the line of sight 
close to the radio core are around 10° iRrichbaum et alJl992h . This 
implies that the precession axis is pointing away from Earth in the 
North; hence 7 = 120° is preferred as the orientation. 

We assume that the jets are precessing anti-clockwise when 
viewed down the precession axis as shown in Fig.0 There are usu¬ 
ally two possible phase angles (<j)) for those bubbles which have 
no constraint on a. The single values have been chosen from the 
distance travelled by the bubble and also from the observed shape; 
for example, the Southern Ghost bubble could have travelled at ei¬ 
ther 70° or 170° to the line of sight. We choose 70° as, if it were 
travelling at 170°, its true age would be much older than calculated 
and the X-ray depression would not be expected to be visible. It 
also appears to be more in cross-section rather than face on. Sim¬ 
ilar arguments have been used to choose the angles for the other 
bubbles. 


6 TIMESCALES 

To calculate the precession timescales we use the change in phase 
between the Ghost and the Outer bubbles, as these are likely to be 
moving rectilinearly and have their ages estimated using the same 
method. In the assumption that the jet is precessing clockwise, then 
for the Northern jet Acj) = 274° and At = 6.1 x 10^ yr; for 
the Southern jet Atf) = 293° and At = 6.2 x 10^ yr. This leads 
to precession timescales of 8.0 and 7.6 x 10^ yr for the Northern 


Tables. Precession Timescales 


North South 


Ghosts 

►Outer 

A(j) = 274° 

Ac/) = 293° 



At = 6.1 

At = 6.2 



'T’prec — 8.01 

'T'prec — 7.62 



Average 

"^prec — 7.8 

Ancient 

—>Ghost 

Ac/) = 61° 

Ac/) = 113° 



At = 3.8 

At = 7.2 



'T'prec — 22.4 

Tpcec = 22.9 



Average ' 

^prec — 22.7 

Outer— » 

Jet 

Ac/) = 201° 

A4> = 249° 



At = 1.5 

At = 1.9 



’^prec — 2.68 

"^prec — 2.75 



Average 

"^prec — 2.7 


Notes: At and Tprec in X10^ yr. 


Table 4. Precession Cycles 




North 

South 


Ghost— 

►Outer 

n = 1 

"^prec — 3.46 

n = 1 

'T'prec — 3.42 


Ancient 

—►Ghost 

n = 1 

Pprec — 3.25 

n = 2 

Pprec =3.11 


Outer— » 

Jet 

n = 0 

’^prec — 2.68 

n = 0 

"^prec — 2.75 


Notes: 

At and Tprec in x 10^ yr. 

n from Tprec = 

360At 

A<^+n360 


and Southern jets respectively, which agree very well. For the other 
precession timescales see Tabled 

Instead of using the Inner bubbles we used those for the VLBI 
jets when comparing to the Outer bubbles. The difference in ages 
between the Inner and Outer bubbles is very short, which would 
give a very rapid precession timescale. In the assumption that the 
VLBI jet is powering the Inner bubbles, it is more sensible to 
choose the parameters for the jets, as the Inner bubbles may still 
be forming. The change in phase angle between the Inner bubbles 
{(j) ~ 178°) and the jets {<j) ~ 182°) is in the correct sense for a 
counter-clockwise precessing jet (in the assumption that the Inner 
bubbles are 4> ^ 178°). 

There is good agreement between the pairs of bubbles on 
the value of the precession timescale, Tprec- There is a large dis¬ 
crepancy between the different bubble pairs. The precession may 
be speeding up, which would explain the decrease in the preces¬ 
sion time; however the acceleration would be very rapid - over 
the course of ~ 1.5 rotations (^^ 600°) the precession timescale 
changes by around a factor of 10, from 0.2 Gyr to 0.02 Gyr. 

It would be expected that the precession timescale would be 
the same over many cycles and so we assume that, as the precession 
timescale for the Outer^Jet is the shortest, then this is the preces¬ 
sion time. The Tprec for the Ghost—»Outer and the Ancient^Ghost 
are much longer, however if an extra precession cycle occurred be¬ 
tween the Ghost and Outer bubbles; i.e. Ai)) = (274+360) = 634° 
for the Northern pair, this would give a precession timescale of 
3.46 X 10^ yr which is much closer to that for the Outer^Iet. For 
the number of cycles required for the other bubbles see Tabled- 

Table 0 shows that, as a result of the large difference in 
timescales for the Ancient Northern and Southern bubbles, an ex- 
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tra precession cycle is required if the precession timescale is fairly 
constant. This could result if we have mis-identified the generation 
the southern Ancient bubble belongs to; i.e. that it is from an even 
older set, and the features of the true Ancient bubble have been lost. 
It may also mean that the bubbles detach at vastly different times; 
the timescales of these old bubbles are calculated using buoyancy 
arguements, so if the northern Ancient bubble detached a long time 
after the southern one, or rose much more slowly, then this may be 
the result. 

The Tprec for the Outer—>Jet are smaller than the other esti¬ 
mates but this can be explained as the bubbles, which correspond to 
the current position of the jets, are still forming and so the preces¬ 
sion timescale for these two may easily be larger. However these 
calculations show that the precession timescale for the SMBH in 
the centre of the Perseus Cluster is around 3.3 x yr. 


7 DISCUSSION 

In the calculation of the precession parameters we have assumed 
that the directions in which the bubbles are observed to be mov¬ 
ing match the directions of the jets at the time when the bubbles 
were just forming. As the timescales for the formation of a bub¬ 
ble (fsound = 1 ^ 2 X 10^ yr) are shorter than the precession 
timescale for the jet (~ 3.3 x 10^ yr), then the bubble will finish 
forming before the jet has moved on (though this would depend on 
the exact formation mechanisms of the bubbles). The travel direc¬ 
tion of the bubbles may therefore indicate the average phase angle 
of the jet during the time of their formation. However, some bub¬ 
bles may rise up in a direction corresponding to the direction of 
the jet towards the end of their inflation phase, and others in the 
direction from the beginning of their formation. There is therefore 
some uncertainty in the phase of the precession corresponding to 
the bubbles and hence in the inferred precession times. 

We have also assumed that the bubbles are not deflected from 
their path once they start to rise up. The bubbles would be expected 
to travel up out of the potential well of the cluster, however a denser 
region of the ICM may deflect the bubbles to different observed 
position angles. 

The “ancient” bubbles that have been used to calculate the 
timescales may not be the remnants of very old bubbles, but just 
chance features in the ICM which have been misinterpreted. How¬ 
ever the fact that the Northern Ancient bubble sits close to the 
low-temperature whorl, and has Ha filaments pointing from the 
centre of the cluster in its direction; the high abundance ridge 
which corresponds to an extension of the radio mini-halo and 
the low-temperature whorl along with th e discovery of some fil¬ 
aments pointing towards the South-West iSanders et alJl200.5ll for 
the Southern Ancient bubble all point towards the fact that these 
features are the remnants of earlier generations of bubbles. 

What can be seen from the phase angles at which the bubbles 
detach is that they do not necessarily detach at the same time. It 
seems likely that the jets produced by the SMBH are of the same 
composition and power, and as such provide the same amount of 
energy at the same rate into the bubbles. From this, in the case that 
the ICM were perfectly uniform, the bubbles would be expected to 
grow at very similar rates and so also detach at very similar times. 
However it is very likely that the ICM is not perfectly uniform, the 
bubbles may expand at slightly different rates, detaching at slightly 
different times. 

The change in the phase of the precession between the detach¬ 
ing of the Ghost and Outer bubbles Is very similar (see Table |3}, 


even with an extra precession cycle. However there is a difference 
of up to 60° between the detaching of the Northern and Southern 
bubbles for a given “generation”, corresponding to around 5 Myr, 
a large fraction of the bubble formation time. In the case of the 
change between the Ancient and Outer bubbles, there is more than 
a complete precession cycle difference between the Southern and 
Northern bubbles (see Table |3- If the Ghost bubbles detached at 
similar times, then the Southern Ancient bubble detached about one 
precession cycle before the Northern Ancient bubble. Whether the 
Southern Ancient bubble belongs to a generation older than that of 
the Northern Ancient bubble, and the true Southern Ancient bubble 
has not yet been identified, or if the bubbles did detach at vastly 
different times is currently unclear. 

This implies that, if the precession is long lived, there are cy¬ 
cles in the precession when no bubbles detach. This may mean that 
the SMBH is dormant and its duty cycle is less than 100 per cent. 
What may also happen is that the jet does not have sufficient power 
to inflate a bubble, but is still supplying energy to the central re¬ 
gions of the cluster by other mechanisms. 

If this were the case, then at times when the jet were less pow¬ 
erful, and possibly less energy were being supplied to the centre the 
ICM may cool and flow inwards for a while. When the jets become 
more powerful, they interact with the cool gas and lift it up out of 
the centre, the result of which we see as the whorl. 

The spirals seen in the optical em ission from NGC1275 
iHoltzman et alll992l : ICarlson et alJll998ll . which are in the same 
sense as the low-temperature whorl are explained if we are seeing 
them from below. The precession of the jet is anti-clockwise when 
seen from above, however the precession axis is tilted away from 
the line of sight, and so the end of a constant length precessing jet 
would appear to move in a clockwise direction, which is the same 
sense as the spirals seen in NGC1275. Their inclination means that 
we see them apparently move clockwise. 

The best fitting parameter s for the precession o f 3C84 are con¬ 
sistent with those presented in iGower et alJ il982h who compared 
models of precessing jets to the observed radio emission. Their 
largest timescale is ~ 10^ yr which is close to the one obtained 
here. The opening angles of the precession have a large range, with 
an average of around ~ 30° which is smaller than the one found 
for 3C84, but not inc onsistent. _ 

Recent work by iLodato & Prinsld j200.5ll show that if black 
hole’s and the accretion disc’s angular momenta were counter- 
aligned to start with, then over the course of what would be one 
precession period the spins would align (in the case where the disc 
angular momentum is large compared to the hole’s). During the 
course of this alignment the spin axis of the black hole traces out 
an arc on the sky. The resultant position angles of the radio jets, 
and hence bubbles, could give the appearance of steady precession 
in the case of the Perseus Cluster. 


8 WHAT COULD CAUSE THE PRECESSION? 

Although a precessing jet model is a possible explanation for the 
observation that the bubbles occur at different position angles in 
the cluster, what is not clear is why the black hole would be pre¬ 
cessing. In a recent comparison between some bipolar Planetary 
Nebulae and the bubbles observed in cl usters, there are two possi - 
ble mechanisms causing the precession iPizzolato & Soke3l200.5ti . 
Firstly there could be a binary black hole at the centre of the 
NGC1275 galaxy, causing the disc of the primary to process iKatzl 
^93)- The second suggestion is that there is an instability in the 
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Table 5. Fluid Dynamics 


Reynolds Number 
Weber Number^ 
Galileo Number 
Strouhal Number 
Froude Number 
Drag Coefficient 
Fluid Parameter 


T^e = 2Uv jv 
We = U'^pr jcF 

G = (\/|po/p-l|9(2r)3)/i/ 

St =fr/{U^) 

Tr = U / 

Co =^^5^ Intel, 

M =g\j4lpa^ 

U((Uoo)) is the (average) bubble velocity, r the bubble radius^, u is 
the kinematic viscosity (/i/p), p and po is the density of the fluid and 
sphere respectively, a is the surface tension, / is the frequency of oscil¬ 
lation of the bubble motion, T^Coo is the asymptotic Reynolds number, 
g is the gravitational accele r ation ^Our definition of We is the squai'e 
of that in iHartunian & Seard ^In some cases an equivalent radius 

re = (3 X Volume/47r)^/^ is used instead. 


disc which warps it and so causes the precession jPrin 

The timescales predicted from these two models are similar; for 

the warped disc model 


^prec 


2 X 10^ yr 


_i M 
“ 109 Afg 


(Equation 4.11 in I Pringld ilQQTt l where a is 
IShakura & Sunva^ L973t viscosity alpha-parameter; 
for the binary model 


the 

and 


/ M \ U 

' a \ 3 / ad N 

-i(l + q)i 

\10^A4qJ '' 


qcosi7 


where q = M 2 /Mi, Mi is the mass of the accreting black hole, 
M 2 is the mass of the object loosing mass, and M = Mi + M 2 
is the total mass of the binary system. The separation between 
the components is a for a circular orbit, ad is the disc radius 
and i9 is the tilt angle betwee n the orbital plane and the disc an¬ 
gle. |pizzoiat^^^ote3 ^200■'5^ prefer the binary-driven precession 
model to the disc instability driven one because there are point 
symmetries in the bubbles in the clusters, assumed to be caused 
by a precessing jet. The warped-disc model would cause a stochas¬ 
tic rather than a regular precession. Both scenarios could produce 
precession times close to those inferred from the bubble directions. 


9 BUBBLE DYNAMICS 

In a still fluid, small gas bubbles rise rectilinearly, however larger 
ones follow zig-zag or spiral paths. This phenomenon has been 
called “Leonardo’s Paradox” as the first recorded observations 
and attempts at explan ations of this behaviour were per formed 
by Leonardo da Vinci iProsnerettil |2004 Ida Vincil| 15151) . Since 
Leonardo’s time there have been many attempts to explain the phe¬ 
nomenon and determine the parameters which cause the change 
between linear and periodic motion. 

The motion of spheres and bubbles in fluids has been studied 
experimentally as well as with numerical simulations. Using solid 
spheres may not at first sight seem appropriate, but if the liquid is 
not pure (“contaminated”) then surfactants on the fluid-bubble in¬ 
terface (partly) immobilise the fluid, similar to the conditions found 
on the surface of a solid object. In this section we review some of 
the work on bubbles rising in fluids to see whether any of these ef¬ 
fects could arise for the bubbles in the Perseus Cluster and cause 
the bubbles to be found at different position angles. In the course 
of these studies the fluids in which the bubbles are moving and the 
motions of the bubbles have been characterised by dimensionless 
quantities which are shown in Tabled 


The fluid in which the bubble moves is described by a dimen¬ 
sionless AI. The behaviour of Cd with TZe depends on the At 
of the fluid. Pluids with 10“^'^ < Ai < 10“® (e.g. filtered wa¬ 
ter, methanol) are “fast” and show the characteristic dip in Cd at 
100 < TZe < 1000, which corresponds to We ~ 1.5. “Slow” flu¬ 
ids (e.g. mineral oil, wate r-glycerin mix) have 10~ ^ < A4 < 10~^ 
and do not show this dip iHartunian & SearJl957^ . 

The critical threshold for the onset of the non-rectilinear mo¬ 
tion depends on the purity of the liquid - the Reynolds number 
being important for contaminated liquids (7?.ecrit ~ 210 ) and the 
Weber number for pure liquids (IVecrit ~ 1.58). 

Numerical simulations show that the wake behind a sphere 
which is held^ed in a uniform flow is axisymmetric up to TZe = 
212. Above this value a planar symmetric wake is found, con¬ 
sisting of two steady counter-rotating threads. At TZe ~ 270 the 
planar symmetric wake become time-dependent; oppositely signed 
stream-wise vortices form a series of loops resembling hairpin 
vortices. Eventually, as TZe continues to increase the wake be¬ 
hind the sphere becomes turbu lent. The exact va lues of TZe de¬ 
pend on the density ratio, po/p Ilennv et all2004i . The results do 
not c hange dramatically i f the sphere is moving freely in a still 
fluid.lde Vries et alH2002h investigate bubble motion as opposed to 
sphere motion and state that in the case of bubbles TZean ~ 600, 
however this discrepancy may aris e from the purity of the fluid. 

Using numerical simulations. Llennv et alJ i2Q04i investigated 
the instabilities in the motions of spheres (rising and falling) in 
Newtonian fluids finding that Gcrit ~ 170, for a density ratio of 
po/p = 0, with a.St = 0.05. They find that a range of behaviours 
occur for bubbles moving with different density ratios and Galileo 
number (see their Lig. 29). Lor a given po/p, the wake behind a 
sphere starts axisymmetric, then as the Galileo number increases 
the motion becomes steady but oblique. Continuing to increase G 
causes oscillations to occur during the oblique ascent, eventually 
causing zig-zag motion. Eventually at C/ ~ 200 the motion be¬ 
comes chaotic. 

There are two types of motion when the wake becomes time 
dependent - helical and zig-zag. The transition between the two 
types of motion i s thought to depend on th e bubble size, observed 
experimentally bv iLunde & Perkin 3 ll 997 h . Larger bubbles deform 
fr om spherical shapes to ellipsoid al ones and numerical simulations 
bv iMougin & Magnaudell L002h show that the transition between 
the two types of periodic motion occurs for an aspect ratio, x = 
2.5. _ 

iHartunian & Sear^ ^19■57^ state that the path instability does 
not occur for small and spherical-cap bubbles, but only for inter¬ 
mediate sized ellipsoidal bubbles, where the motion ranges from 
oscillations to gentle rocking about the direction of motion. It is 
possible for the zi g-zag motions to cause d eflections from a verti¬ 
cal path by to 50° lUunde & Perkin Jl9vili which is similar to the 
position angles of the North-Western Ghost and the South-Western 
Ancient bubbles. 


9.1 Spherical Cap Bubbles 

IWegener & Parlan^^l973^ show that Ar = 1.00±0.05 for spher¬ 
ical cap bubbles. In the assumption that this relation continues at 
much larger volumes, then for the North-wester n Ghost bubble, the 
buoyancy velocity calculated during the work of iDunn et alJ l2003i 
is ~ 3.6 X 10^ cm and the radius of the bubble is ~ 14kpc. 
Using the enclosed mass and i7dist to obtain an estimate on the 
gravitational field, q, gives t Fr ~ 0.92. Therefore, if the relation in 
IWegener & Parians il973^ holds up to these bubble sizes, then it 
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is likely that the observed X-ray decrement is, and behaves like, a 
spherical cap bubble. 

In this case it could be a s sumed to follow the behaviour 
observed by iHartunian & Sear^ ^9^), i.e. it rises rectilinearly. 
Therefore it could be assumed that the other Ghost and Ancient 
bubbles would have travelled linearly out from the cluster after they 
had attained the spherical-cap shape. 


9.2 Cluster Bubbles 

The bu bbles in the Perseus Cluster as analysed by iFabian et alJ 
j2Q03dl inferred that the ICM was viscous, and as such the 
Reynolds number would be expected to be less than 1000. The 
kinematic viscosity inferred from this Reynolds number and the 
North-western Ghost bubble is 4 x 10^^ cm^ s“^. Therefore it 
could reasonably be assumed that the ICM is a “slow, contami¬ 
nated” fluid and as such the Reynolds number rather than the We¬ 
ber number would determine the onset of path instability in the 
bub bles. _ 

IWegener & Parian^ ^1973^ state that spherical cap bubbles 
rise steadily with no periodic motion, and so the path instability 
would apply to the newly formed bubbles as they would be spher¬ 
ical or ellipsoidal. The critical Reynolds number for the onset of 
the path instability is ~ 200. Therefore, as the upper limit on the 
Reynolds number of the ICM is ~ 1000, the bubbles could have 
been oscillating about their direction of motion shortly after their 
creation. As the bubbles’ shape evolved into a spherical cap their 
oscillations reduced and they may then end up travelling at a dif¬ 
ferent position angle than that of the jet. 

For the South-western Outer bubble, using the kinematic vis¬ 
cosity from iFabian et alJ i2003dl . and assuming that the density 
ratio po/p ~ 0, then the Galileo number is ~ 1140, which is in the 
regime of chaotic motion. This is unlikely as there are indicators in 
the X-ray and Ha emission that the motion of the bubbles has been 
fairly steady. As such the Galileo number would be expected to be 
much lower, < 200. This implies that the Reynolds number, which 
has been used to estimate the viscosity, would be around 200, giv¬ 
ing a viscosity of 2 x 10^® cm^ s“^. There is no large change in 
these values if po/p < 1 were used instead. 

The mechanism which causes the onset of path instability in 
the rising bubbles is currently uncertain. However, if the onset of 
the path instability occurs at the same values of the Galileo and 
Reynolds numbers at the scales observed in clusters of galaxies, 
then, as the calculated values for these numbers fall in the range 
where the bubble motion would be expected to be non-rectilinear, 
the bubbles could oscillate as they rise. This motion, however, 
would not last for long as the bubbles appear to evolve into spheri¬ 
cal caps (which show no oscillatory motion, see Section lgTl com¬ 
paratively quickly, before complete oscillations would have taken 
place. The position angles of the inner bubbles are such that the de¬ 
viations from a linear rise would have to start very quickly, almost 
immediately after the bubbles detached. 

The results of the experiments discussed in the above sections 
have all been performed on the Earth where the gravitational po¬ 
tential can be regarded as 1-dimensional. In the environment of the 
cluster the gravitational potential is spherically symmetric. What 
effect this would have on bubbles which are oscillating about their 
direction of motion is not clear, and probably depends on rela¬ 
tive sizes of the bubble to the cluster potential. Bubbles moving 
obliquely in experiments on Earth may move in spirals in a cluster, 
which could be an explanation for the low temperature structure 


observed in the temperature maps of ISanders et alJ ^2Q04^ if the 
bubbles did not evolve into a spherical-cap shape. 

However this does depend on there being sufficient time for 
the oscillations to occur before the bubble evolves into a spherical 
cap shape; the Southern Outer bubbles’ shape as seen in the X-rays 
already appears similar to that of a spherical cap. If the Strouhal 
number relation also holds up to these bubble volumes, then / « 
10“^® ~ 3 Gyr“^. This implies that there would be no chance 

for any oscillation before the bubble shape evolves into a spherical 
cap. 

The radial arrangement of the Ha filaments in the cluster 
strongly suggest that the bubbles move in a linear fashion, as os¬ 
cillating bubbles, or bubbles moving in spirals would be expected 
to drag the filaments with them as they rise. 


10 CONCLUSIONS 

The observations of the bubbles in the Perseus Cluster show that 
they are at different position angles with respect to the cluster cen¬ 
tre. We have outlined several interpretations for this behaviour, in¬ 
cluding a jet which is precessing or moving around. Eor the pre- 
cessing jet we find that the best fit precession parameters are that 
the projected precession angle offset is 10° westwards of North; 
the precession axis is tilted 120° away from the line of sight and 
the precession opening angle is 50°. Using the timescales for the 
bubbles the precession timescale is around 3.3 x 10^ yr. We review 
some work on the motion on bubbles rising in fluids and the onset 
of path instability. The conditions of the motion in the cluster are 
such that the oscillations about the direction of motion are possi¬ 
ble for the young bubbles, but the timescales are of order 1 Gyr, 
by which time the bubbles have evolved into spherical caps, which 
do not exhibit this type of motion. A steady precession or motion 
of the jet may be responsible for the whorl seen in the temperature 
structure. Alternatively, we note that the whorl may be an artifact of 
a partial inner circle and a partial outer ellipse, connected together 
in the SE. The structures at the centre of the Perseus cluster provide 
an enigmatic view of the past history of this region. 
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